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Abstract.  A  new  computational  method  for  characterizing  the  relationship  between  surface  crystallography  and 
electrical  conductivity  in  anisotropic  materials  with  low  frequency  electromagnetic  techniques  is  presented.  The  method 
is  discussed  from  the  standpoint  of  characterizing  the  orientation  of  a  single  grain,  as  well  as  characterizing  statistical 
information  about  grain  ensembles  in  the  microstructure.  Large-area  electron  backscatter  diffraction  (EBSD)  data  was 
obtained  and  used  in  conjunction  with  a  synthetic  aperture  approach  to  simulate  the  eddy  current  response  of  beta 
annealed  Ti-6A1-4V.  Experimental  eddy  current  results  are  compared  to  the  computed  eddy  current  approximations 
based  on  electron  backscatter  diffraction  (EBSD)  data,  demonstrating  good  agreement.  The  detectability  of  notches  in 
the  presence  of  noise  from  micro  structure  is  analyzed  with  the  described  simulation  method  and  advantages  and 
limitations  of  this  method  are  discussed  relative  to  other  NDE  techniques  for  such  analysis. 

Keywords:  eddy  current,  micro  structure,  hexagonal,  anisotropic 
PACS:  81.70.-q,  81.05.Xj 


INTRODUCTION 

One  of  the  main  confounding  factors  in  using  the  methods  of  eddy  current  testing  (ECT)  to  detect  small  (<  0.1  mm) 
cracks  in  engineering  anisotropic  alloys  is  the  micro  structure  of  the  material.  When  the  change  of  impedance  due  to 
the  presence  of  the  crack  is  on  the  same  order  as  the  change  due  to  perturbations  in  crystallographic  orientation,  the 
crack  signal  gets  lost  in  what  is  often  perceived  as  the  surrounding  “grain  noise”  [1].  This  situation  is  illustrated  in 
Fig  1.  The  same  sized  electrical  discharge  machined  (EDM)  notch  was  scanned  in  two  different  materials.  The 
main  difference  between  the  two  images  is  that  the  material  in  the  image  on  the  left  is  isotropic  vs.  the  anisotropic 
material  on  the  right.  Several  methods  could  be  used  to  enhance  the  signal  from  the  notch  in  FIGURE  1(b).  Any  of 
several  image  processing  routines  could  be  used  to  identify  and  subtract  the  features  caused  by  the  micro  structure 
from  that  due  to  the  presence  of  damage.  This  would  offer  some  sensitivity  to  the  notch  or  it  could  provide  a 
method  to  reduce  the  noise.  This  would  be  a  similar  approach  to  using  spatial  Fourier  transforms  to  extract  features 
from  images.  These  methods  do  offer  some  enhanced  sensitivity  to  the  notch  signal,  but  in  the  process  they  make 
parameter  estimation  much  more  involved.  For  NDE  purposes,  the  ideal  situation  is  one  in  which  this  data  can  be 
used  in  its  raw  form  to  estimate  characteristics  of  the  damage.  For  this  type  of  analysis,  a  thorough  understanding  of 


(a)  (b) 

FIGURE  1.  Eddy  current  scans  taken  over  the  same  sized  crack  in  (a)  AL  2024  and  in  (b)  Ti-6A1-4V  (P  annealed)  [1] 
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FIGURE  2.  Image  of  the  sample  analyzed  in  this  work.  The  grains  are  clearly  large  enough  that  the  coil  can  resolve  them. 

the  physics  involved  is  extremely  helpful.  To  this  end,  this  work  is  mainly  concerned  with  proof  of  the  theory  that 
this  noise  is  related  to  microstructure  and  in  providing  a  potential  forward  model  approximation  for  this  type  of 
problem. 


EXPERIMENTS 

Experiments  were  performed  initially  to  determine  the  source  of  the  material  noise  and  then  to  provide  sample 
information  to  the  forward  models  for  verification.  An  appropriate  sample  was  chosen  for  the  purposes  of  the  work, 
and  then  several  different  experimental  techniques  were  used  to  collect  information  about  the  microstructure.  The 
different  experimental  data  sets  were  compared  for  verification. 

Sample 

In  order  to  illustrate  and  exacerbate  the  role  of  grain  orientation  on  conductivity,  thereby  demonstrating  the 
source  of  microstructure-induced  noise,  a  coarse-grained  beta-annealed  Ti-6A1-4V  sample  was  selected  for  this 
work.  The  grains  were  sufficiently  large  that  the  diameter  of  a  standard,  commercially  available  ECT  probe  was  less 
than  the  typical  grain  diameter  and  therefore  the  signal  at  a  given  could  be  due  to  a  single  grain.  It  is  noteworthy 
that  the  large  grain  dimensions,  up  to  -1  cm,  in  this  sample  imply  a  large  radius  of  curvature  of  the  colony 
boundaries,  much  larger,  in  fact,  than  the  depth  contributing  to  the  eddy  current  signal  and  thus  complications  due  to 
subsurface  variations  in  orientation  that  may  be  present  in  fine  grained  materials  were  not  a  source  of  variability 
here.  The  optical  images  of  the  sample  shown  in  Fig.  2  illustrate  this  point  more  clearly.  Ti-6A1-4V  consists 
primarily  of  the  hexagonal  close  packed  a-phase  (-93%  volume  fraction)  with  the  remainder  balanced  by  the  body 
center  cubic  p-phase.  The  a-phase  exhibits  isotropy  in  the  basal  plane,  but  most  physical  properties  are  moderately 
anisotropic  and  depend  on  angle  relative  to  the  [0001]  (c-axis)  direction.  The  variation  in  conductivity  from  the 
basal  plane  to  the  c-axis  is  nominally  about  6%  difference,  and  the  conductivity  in  the  basal  plane  is  assumed  to  be 
on  the  order  of  2.4xl06  S/m  [2]. 


Electron  Backscatter  Diffraction 


Since  most  properties  vary  as  a  function  of  orientation,  and  the  cost  of  these  measurements  is  becoming  more 
affordable,  large-area  microstructure  characterization  involving  measurement  of  local  orientation  via  EBSD  is 
becoming  increasingly  commonplace  (e.g.  [3]).  Crystallographic  orientation  in  an  orthonormal  basis  can  be  fully 


FIGURE  3.  (a)  The  hexagonal  crystal  in  an  orientation  in  line  with  the  sample  axes,  (b)  rotated  an  angle  (p±  about  the  initial 
z-axis,  (c)  rotated  an  angle  6  about  the  new  x-axis,  and  (d)  rotated  an  angle  cp2  about  the  newest  z-axis. 
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FIGURE  4.  Crystal  orientation  map  of  the  coarse-grained  Ti-6A1-4V  sample.  Here,  color  corresponds  to  crystallographic 
orientation.  The  large  grain  in  the  center  is  outlined  for  reference. 

described  by  a  set  of  Euler  angles  ((Pi,0,(p2)  in  Bunge’s  active  form,  as  shown  in  Fig.  3.  The  Euler  angles 
essentially  define  the  rotations  necessary  to  bring  the  sample  axes  into  coincidence  with  an  arbitrarily  oriented 
crystal  [Error!  Reference  source  not  found.].  In  this  work,  a  direct  comparison  between  eddy  current  images  and 
the  actual  sample  cannot  be  made  without  this  orientation  information.  Orientation  data  were  collected  with  electron 
backscatter  diffraction  using  a  Phillips  XL30  scanning  electron  microscope  (SEM)  operated  at  an  accelerating 
voltage  of  20  kV  at  a  probe  current  of  50  nA.  The  resulting  data  is  shown  in  Fig  4. 

Eddy  Current  Experiments 

A  commercially  available,  2MHz  center  frequency  probe  was  used  to  collect  ECT  data  from  the  large  grained 
sample  shown  above.  The  probe  was  connected  to  a  Nortec  2000D  eddy  current  scope  and  raster  scanned  over  the 
entire  specimen.  This  combination  of  instrument  and  probe  has  been  shown  to  produce  detectable  signals  from 
conductivity  changes  of  below  6%,  which  is  the  maximum  change  that  would  be  expected  from  the  sample  (i.e.  the 
percent  difference  between  the  c-axis  and  the  basal  plane  conductivities).  To  achieve  such  sensitive  measurements, 
one  crucial  step  is  to  determine  the  probe’s  sensitivity  to  conductivity  vs.  liftoff.  A  calibration  curve  for  this  probe 
generated  with  physical  models  [5]  can  be  seen  in  Fig.  5.  The  conductivity  range  that  the  sample  is  in  is  close  to 
1.75  %IACS,  which,  upon  inspection  of  the  curve  in  that  range,  means  that  liftoff  is  not  orthogonal  to  conductivity 
changes.  However,  with  conductivity  changes  on  the  order  of  6%  [2],  the  magnitude  change  is  close  enough  to 
linear  to  use  it  as  the  impedance  measurement.  Liftoff  was  controlled  by  using  the  probe  in  contact  with  the 
specimen.  An  image  of  the  magnitude  of  impedance  above  the  sample  is  shown  in  Fig.  6.  This  image  shows  some 
of  the  grains  that  can  be  seen  in  the  OIM  image  indicating  that  the  ECT  method  is  capable  of  distinguishing  the 
grains,  and  that  the  change  in  the  impedance  of  the  coil  is  due  mainly  to  the  conductivity  change  from  one  region  to 
another. 


APPROXIMATION  METHOD 

Knowing  that  this  variation  in  impedance  measurements  over  Ti  and  its  alloys  is  a  conductivity  variation  due  to 
microstructure,  a  quick  and  accurate  approximation  method  for  forward  models  of  this  phenomenon  can  be  obtained. 
This  is  useful  if  statistical  quantification  of  the  micro  structure  is  to  be  accomplished  with  ECT  methods.  The  key  to 
the  approximation  lies  in  several  assumptions  that  will  be  enumerated  here  to  clarify  the  method  as  well  as  provide 
potential  avenues  to  refine  the  approximation. 


Impedance  Plane  Plot  for  2MHz  Coil 


FIGURE  5.  Calibration  curve  for  the  eddy  current  probe.  The  red  line  represents  change  in  conductivity,  and  the  blue 

represents  a  change  in  liftoff  at  a  fixed  conductivity. 
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FIGURE  6.  ECT  scan  of  the  sample 

Assumption  1:  Uniform  Electric  Field 

This  assumption  states  that  underneath  the  coil,  the  electric  field  is  uniform  inside  the  conductor.  This  is  an 
assumption  since,  as  can  be  seen  in  any  mathematical  analysis  of  eddy  current  problems  [6,  7],  the  electric  field  in 
the  conductor  is  dependent  on  the  conductivity  of  the  sample.  However,  this  assumption  hinges  on  the  fact  that  the 
conductivity  changes  are  so  small  that  the  variation  of  electric  field  from  one  grain  to  the  other  is  small.  The 
definition  of  “small  enough”  here  is  not  analyzed,  and  in  fact  this  assumption  might  be  fairly  stringent.  However, 
since  the  goal  is  to  produce  an  approximation  of  an  eddy  current  scan,  this  assumption  will  be  enforced. 

Assumption  2:  No  Changes  in  Depth 

The  second  assumption  is  that  the  grains  that  we  see  in  the  EBSD  map  are  “pillars”,  or  that  they  are  constant  in 
the  depth  of  the  sample.  This  is  clearly  a  strict  and  rather  limiting  assumption.  The  skin  depth  of  the  eddy  currents 
in  this  material  at  2MHz  is  around  230pm.  Common  commercially  used  Ti  alloys  have  micro  structures  with  grain 
sizes  in  the  range  of  10-20pm  meaning  that  the  present  analysis  would  be  confounded  by  non  surface-connected 
grains  which  would  contribute  to  the  eddy  current  response.  In  this  case,  however,  the  large  grain  size  and 
commensurate  large  grain  boundary  curvature  implies  that  the  region  of  constant  orientation  far  exceeds  the  depth  of 
the  eddy  currents.  Optical  inspection  of  both  sides  of  the  sample  indicated  that  many  of  the  grains  actually  spanned 
the  entire  thickness  (~5  mm)  validating  that  assumption  that  the  EBSD  measurements  are  representative  of  the 
orientations  at  least. 


Details  of  the  Method 

With  the  assumptions  clear,  the  approximation  routine  can  now  be  elaborated  upon.  The  first  thing  that  must  be 
understood  is  that  the  conductivity  in  the  material  is  different  for  electric  fields  of  different  direction.  The 
constitutive  relation  for  induced  current  density  in  the  presence  of  an  electric  field  is  [8]: 

J=  gE  (1) 

J  is  the  current  density,  o  is  the  conductivity  tensor,  and  E  is  the  electric  field.  If  we  assume  an  electric  field  in  the 
direction  of  a  unit  direction  vector,  r,  and  assign  it  some  arbitrary  magnitude,  say  Emag ,  then  with  some  steps  of 
linear  algebra,  the  magnitude  of  the  current  density  induced  can  be  expressed  as: 

Jmag  =  Emag(rT  (J2^1/2  (2) 

This  expresses  the  magnitude  of  the  current  density  as  a  function  of  an  electric  field  in  an  arbitrary  direction  and  the 
conductivity  tensor  of  the  sample.  The  current  density  here  is  not  necessarily  in  the  direction  of  the  electric  field, 
but  its  magnitude  is  correct.  We  can  define  an  effective  conductivity  as  oeff  =  (rrcr2r)1/2.  When  the  electric  field 
is  in  the  direction  of  one  of  the  principal  axes  of  the  conductivity  tensor,  the  conductivity  comes  out  as  the  single 
crystal  conductivity  in  that  direction,  as  expected.  Now,  for  an  arbitrary  set  of  rotations  given  Euler  angles,  a 
rotation  matrix  for  the  conductivity  tensor  can  be  defined  as  [4] : 
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T 


(3) 


R  = 


clc3  ~  C2S1S3 

C3S 1  +  c1c2s3 

S2S3 


~cls3  ~  C2C3S1 
C1C2C3  ~  sls3 
C3S2 


sls2 

-ClS2 

c2 


C  represents  cosine  and  s  represents  sine,  and  the  subscripts  1,  2,  and  3  represent  the  operations  on  Euler  angles  <p1? 
0 ,  and  cp2  respectively.  Given  a  new  direction,  Rr ,  for  the  electric  field,  the  electric  field  can  be  expressed  as: 

E  —  EmagRr  (4) 

Again,  we  write  the  constitutive  relation  for  the  newly  aligned  region  as: 

J  —  &E  —  EmagoRr  (5) 

We  can  carry  out  the  same  operations  for  an  electric  field  in  an  arbitrary  direction  and  arrive  at  a  final  expression  for 
the  magnitude  of  the  current  density  as: 


Jmug  =  ElnafJ(rTRTa2Rr)1/2  (6) 

where  we’ve  used  the  fact  that  R_1  =  RT .  The  effective  conductivity  in  this  expression: 

aejrjr  =  (rT  Rt a2  Rr)1/2  (7) 

represents  the  conductivity  in  the  direction  of  the  electric  field. 

With  this  definition  of  effective  conductivity,  an  expression  for  the  conductivity  of  a  point  below  the  coil  can  be 
derived.  The  first  assumption  implies  that  the  field  below  the  coil  and  in  the  sample  flows  tangential  to  the  coil. 
This  is  a  similar  assumption  to  that  made  when  solving  an  axisymmetric  problem  with  scalar  conductivities  as  in  [5]. 
With  that  in  mind,  the  direction  of  the  electric  field  at  a  point,  below  the  coil  can  be  calculated  by  drawing  a  line 
to  the  point  and  then  drawing  the  normal  line  to  that  point  in  the  direction  opposite  the  flow  of  currents  in  the  coil 
(due  to  Lentz’s  law,  though  the  direction  doesn’t  matter  in  this  application).  The  schematic  of  this  is  shown  in  Fig. 
7.  Setting  the  coil  center  as  the  origin  (0,0)  and  calling  the  vector  to  the  point  the  equation  for  the  tangential 
direction,  rh  can  be  written  as: 


[0  0  1]  x  dt  _  [~dy  dx  0] 
\d[\  "  \d[\ 


(8) 


This  value  of  rh  along  with  the  Euler  angles  at  the  point  are  then  use  in  equation  7  to  calculate  the  effective 
conductivity  at  that  point.  Since  the  OIM  data  is  a  regular  grid  of  points,  averaging  every  point  below  the  coil  gives 
what  is  similar  to  an  area- weighted  average  of  orientation-dependent  conductivities.  This  average,  when  plotted  for 
every  coil  position,  gives  an  estimate  of  what  an  ECT  scan  would  look  like  above  an  anisotropic  material.  This 
essentially  amounts  to  a  synthetic  aperture  approach  for  eddy  current  scans  above  a  conductive  medium  with  subtle 
variations  in  conductivity. 


FIGURE  7.  The  direction  of  the  electric  field  used  in  effective  conductivity  calculations  at  a  point,  pi,  below  the  coil. 
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(b) 


FIGURE  8.  Images  of  (a)  an  actual  ECT  scan  of  the  sample  and  (b)  a  simulated  ECT  scan. 

COMPARISONS  AND  VERIFICATION 

To  test  the  synthetic  aperture  routine,  the  Euler  angles  that  were  collected  from  the  sample  were  fed  into  the 
numerical  method  and  a  simulated  eddy  current  scan  was  performed.  The  images  are  shown  back-to-back  in  Fig  8. 
There  are  some  features  that  are  clearly  visible  in  both  images,  but  other  regions  are  very  difficult  to  compare.  The 
cause  of  this  is  partly  the  fact  that  the  first  assumption  implies  that  the  results  of  the  simulation  are  not  physically 
accurate,  and  that  the  3D  grain  morphology  is  missing  from  the  simulation.  Again,  since  the  grains  are  mostly  large 
enough,  the  second  assumption  is  not  as  crucial,  but  it  could  be  relevant  if  the  angle  of  the  grain  boundary  is  drastic 
in  depth.  Further  validation  to  estimate  the  effect  of  the  assumptions  would  require  3D  sectioning  experiments  and 
the  use  of  a  full  numerical  solution  such  as  finite  element  method  (FEM),  boundary  element  method  (BEM),  or 
volume  integral  method  (VIM).  As  of  writing  this  paper,  these  experiments  have  not  been  performed. 

On  the  other  hand,  this  type  of  verification  has  been  demonstrated  to  a  lesser  extent  with  FEM  methods.  A 
problem  with  the  geometry  and  conductivities  shown  in  Fig.  9  was  analyzed  with  COMSOF  Multiphysics.  The 
details  of  this  model  can  be  found  in  [9].  The  two  different  conductivity  regions  were  simulated  using  a  different 
rotation  matrix  for  each.  The  FEM  simulation  gives  changes  in  impedance  of  the  coil,  whereas  the  synthetic 
aperture  routine  give  changes  in  conductivity.  However,  it  was  already  shown  that  this  relationship  can  be 
considered  linear  with  such  small  changes  in  conductivity,  so  the  results  from  each  were  normalized  and  plotted  in 
Fig.  10.  The  synthetic  aperture  routine  shows  very  good  agreement  with  the  FEM  results. 

CONCLUSIONS,  DISCUSSION,  AND  FUTURE  WORK 

From  the  data  provided  in  this  paper,  it  is  reasonable  to  conclude  that  noise  in  the  scans  shown  is  likely  due  to 
microstructure  changes  in  an  anisotropic  material.  This  conclusion  has  been  reach  by  other  authors,  but  no  direct 
comparison  between  OIM  and  eddy  current  scans  has  been  made.  The  OIM  maps  of  these  scans  offer  insight  into 
the  actual  composition  of  the  sample,  as  well  as  give  quantitative  information  that  can  be  used  in  the  context  of 
forward  models.  The  approximation  method  shown  here  has  potential  to  expedite  these  forward  model  solutions  by 
providing  the  ECT  information  over  the  micro  structure.  This  bodes  well  for  future  inversion  schemes  that  may  take 
advantage  of  the  physical  intuition  of  anisotropic  material  noise.  Incorporation  of  this  into  parameter  estimation 
routines  remains  as  future  work.  First  and  foremost,  a  very  thorough  investigation  of  the  validity  of  the  assumptions 
made  must  be  undertaken  and  the  approximation  method  must  be  validated  with  real  experimental  data. 


Coil 

n 

(T1  =  Ri& 

<J2  =  ^2° 

FIGURE  9.  Image  of  the  geometry  analyzed  with  the  FEM  model  as  well  as  the  approximation  method. 


6 

Approved  for  public  release;  distribution  unlimited. 


FIGURE  10.  Results  from  the  FEM  model  are  shown  in  red  dots,  and  results  from  the  synthetic  aperture  routine  are  shown 

as  a  blue  line 
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